pub fn new_birthday_probability(n: u32) -> f64 {
    1. - (1..n).map(|x| (1. - (x as f64 / 365.))).product::<f64>()
}
